Comparison of numerical solution strategies for gravity field recovery from GOCE SGG observations implemented on a parallel platform

The recovery of a full set of gravity field parameters from satellite gravity gradiometry (SGG) is a huge numerical and computational task. In practice, parallel computing has to be applied to estimate the more than 90 000 harmonic coefficients parameterizing the Earth’s gravity field up to a maximum spherical harmonic degree of 300. Three independent solution strategies, i.e. two iterative methods (preconditioned conjugate gradient method, semi-analytic approach) and a strict solver (Distributed Non-approximative Adjustment), which are operational on a parallel platform (‘Graz Beowulf Cluster’), are assessed and compared both theoretically and on the basis of a realistic-as-possible numerical simulation, regarding the accuracy of the results, as well as the computational effort. Special concern is given to the correct treatment of the coloured noise characteristics of the gradiometer. The numerical simulations show that there are no significant discrepancies among the solutions of the three methods. The newly proposed Distributed Nonapproximative Adjustment approach, which is the only one of the three methods that solves the inverse problem in a strict sense, also turns out to be a feasible method for practical applications.